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Abstract 

The azimuthal asymmetry of minijet system produced at the early stage of 
nucleon-nucleon and nuclear collisions in a central rapidity window is studied. We 
show that in pp collisions the minijet transverse energy production in a central 
rapidity window is essentially unbalanced in azimuth due to asymmetric contri- 
butions in which only one minijet hits the acceptance window. We further study 
the angular pattern of transverse energy flow generated by semihard degrees of 
freedom at the early stage of high energy nuclear collisions and its dependence 
on the number of semihard collisions in the models both including and neglecting 
soft contributions to the inelastic cross section at RHIC and LHC energies as well 
as on the choice of the infrared cutoff. 
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1 Introduction 



Minijet physics is one of the most promising apphcations of perturbative QCD 
to the analysis of processes with multiparticle production. It addresses a crucial 
question of how many (semi)hard degrees of freedom can be available in a given 
event. The approach is based on the idea that some portion of transverse energy 
is produced in the semihard form, i.e. is perturbatively calculable because of 
the relatively large transverse momenta involved in the scattering, but, due to 
parametrically strong hadronization effects, can not observed in the form of cus- 
tomary well coUimated hard jets distinctly separable from the soft background 
This mechanism operates at the early stage of the collision and, when relevant, 
determines the characteristics of the primordial transverse energy flow. 

Creation of many (semi)hard degrees of freedom corresponds to a new physi- 
cal situation characterized by nontrivial, possibly kinetic or even hydrodynamic, 
phenomena occurring at the parton level at the early stages of high energy colli- 
sion. Of special relevance are here the ultrarelativistic heavy ion collisions, where 
one would expect that creating a dense system of (semi) hard degrees of free- 
dom in the volume much larger than, e.g., the proton one, possible, thus making 
the application of concepts borrowed from macroscopic physics natural. Recent 
critical discussion of this field can be found in 

Minijet physics is an actively developing field. Reviews on the subject con- 
taining a large number of references are, e.g., 0, and [Q. Several approaches 
have been considered with an aim of providing a quantitative description of pri- 
mordial parton system produced at the earliest stage of, e.g., high energy heavy 
ion collisions. The conceptually simplest one is based on the standard formal- 
ism of coUinearly f actor ized QCD at small parton densities, see 0, 0, 0, 
@. Here one operates with a single hard parton-parton scattering in a given 
hadron-hadron collision, so that standard QCD structure functions can be used 
in computing the probability of generating a pair of partons with certain kine- 
matical characteristics. 

This approach has a natural generalization, in which multiple binary parton- 
parton collisions in the given hadron-hadron one are considered, provided some 
ad hoc distribution in the number of these collisions is chosen, see e.g. [IC]. 
This also allows to construct geometrically motivated scheme of unitarizating 
the semihard contribution to the inelastic cross section of hadron scattering as 
described, e.g., in 0. 



Starting from |jTT[, nonlinear QCD effects in relation with describing the early 
stages of heavy ion collisions drew progressively more and more attention. New 
results were obtained within the approach to minijet production based on the 
quasiclassical treatment of nuclear gluon distributions within a framework of 
McLerran-Venugopalan model []T2|, see [|r^, [jT3|, |TB| and |T^. First non- 



perturbative results on gluon production are now available |jT8| , see also and 
Recently a nonperturbative model for gluon production in heavy ion colli- 
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sions based on the physical concepts of McLerran-Venugopalan model and a cor- 
responding kinetic equation describing the evolution of primordial gluon system 
were discussed in [^. For pedagogical introduction to this rapidly developing 



field see the lectures |^ and p3| . 

A notable feature of the physical phenomena related to the collective behavior 
of multiparton systems is their genuine event-by-event nature, so that many usual 
tools used in high energy physics, such as inclusive distributions, are becoming 
less helpful. Thus the analysis of event-by-event variations of the quantities 



sensitive to the collective dynamics is very important, see. e.g. p4[, |2^, [ P6[] , 
and references therein. 

The description of the primordial parton configuration should provide infor- 
mation on the event-by-event pattern of the parton system, in particular on the 
number of perturbative parton-producing interactions, which, to a large extent, 
determines the initial parton density and other kinematical characteristics. In 
particular, the discrete nature of parton production in phase space, as described 
by finite order QCD calculations, can give rise to primordial event-by-event angu- 
lar asymmetries of the parton flow. The fate of the primordial angular asymme- 
tries depends on the relevant dynamics (is the evolution of the produced system 
long enough to wash them out, can they be frozen and directly relevant to the 
observed hadronic flow, etc). In any case the first problem to look at is to study 
the primordial parton system before the reinteraction of partons sets in. 

The aim of this paper is to study the characteristics of the initial minijet- 
induced transverse energy flow in nucleon-nucleon and nucleus-nucleus collisions 
within the framework of minijet production scenario based on collinearly factor- 
ized QCD, p and 0. In particular, we shall analyze the fiuctuational azimuthal 
imbalance in the minijet transverse energy flow due to a discrete nature of trans- 
verse energy production through basic QCD hard scattering. We can expect that 
the effect will be essentially sensitive to the number of semihard scatterings. In 
what follows we shall see that this is indeed the case. 

Below we study the event-by-event inhomogenities in the azimuthal distribu- 
tion of minijets following from the basic asymmetry of minijet transverse energy 
production into a finite rapidity window in pp collisions. The nuclear collisions 
are described by a geometric model in which they are considered as a super- 
position of basic nucleon-nucleon ones. The azimuthal asymmetry of the minijet 
system will be specified in terms of (transverse) momenta only, as calculated in 
the conventional S-matrix field theory formalism without referring to the coor- 
dinates of partons and making no assumptions on the structure of the contri- 
butions of higher order in QCD coupling constant in describing the transverse 
energy production in the elementary nucleon-nucleon collision. This analysis can 
be extended to the next-to- leading order (e.g. along the lines of [|28[) due to 
the infrared stability of the considered distributions which are of energy-energy 
correlation type. 

The analysis of the event-by-event pattern of the initial minijet generated 
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transverse energy flow was first presented in , where a HIJING generated 
list of partons with specified coordinates and momenta was used to compute a 
coarse - grained energy density and velocity field at RHIC energy y/s = 200 GeV. 
The resulting distributions turned out to be highly irregular and similar to the 
ones occurring in turbulent flows. Note that besides the parton-parton scattering 
described by collinearly factorized QCD considered in the present paper, HIJING 
makes explicit assumptions on the structure of higher order contributions (uni- 
tarization), contribution from initial and final state radiation and structure of 
parton system in coordinate space (thus going beyond standard S-matrix formal- 
ism). The existence of asymmetry due to imbalanced particle production from 
minijets into a finite acceptance was mentioned in 0. 
The outline of the paper is as follows. 

In the second section we analyse the basic mechanism for producing az- 
imuthally symmetric and asymmetric configurations in the restricted phase space 
domain, which in the considered case will be a unit central rapidity window, in pp 
collisions. We calculate in the leading twist (lowest order in parton density) and 
leading order collinear factorization scheme the relative weights for symmetric 
(two-jet) and asymmetric (one-jet) contributions to the transverse energy pro- 
duction cross section for RHIC l^/s = 200 GeV) and LHC (v^ = 5500 GeV) for 
underlying nucleon-nucleon collisions. 

In the third section the computed contributions to azimuthally symmetric and 
asymmetric components of the pp minijet transverse energy production into a unit 
central rapidity window are used in calculating the asymmetry of transverse en- 
ergy production in heavy ion collisions, where the nuclear collision is described as 
a superposition of the nucleon-nucleon ones in a geometrical approach of 0. We 
study the azimuthal asymmetry for RHIC and LHC energies for central collisions 
for two dynamical scenarios. In the first scenario the transverse energy produc- 
tion is assumed to occur through two physically different mechanisms, the soft 
one and the (semi)hard one. As our aim is to study the transverse energy flow at 
early collisions stage related to the semihard degrees of freedom, the contribution 
of soft interactions will be accounted for only in determining the relative yield 
of semihard contribution. In the second scenario, which can become realistic at 
LHC energies, all primordial transverse energy production is assumed to occur 
through the semihard mechanism. 

In the last section we discuss the results and formulate the conclusions. 



2 Azimuthal pattern of minijet production in pp 
collisions 

The mechanism responsible for transverse energy production in the leading order 
in QCD perturbation theory is an elastic two-to-two parton-parton scattering. 



4 



Its cross section is given by the standard coUinearly factorized expression 

Xifixi,pl)-^X2fix2,pi), (1) 



dp\dyidy2 dp] 

where xf{x,p\) is a parton structure function, xi^2 = p±{e^^^ + e^^^)/\^ are 
the fractional longitudinal momenta of the produced partons and da/dp\ is a 
differential cross section of elastic parton-parton scattering. In the following we 
will be specifically interested in the transverse energy production into some given 
(central) rapidity interval ymin < Z/1,1/2 < Vmax- Operationally the transverse 
energy E± deposited in this window by the two scattered partons is defined as [] 

E± = Pl 6{ymin <yi< Vmax) + P2 OiVmin < Z/2 < Vmax) (2) 

In the following we shall confine ourselves to considering the central rapidity 
interval ymin = —0.5 < y < ymax = 0.5 and stay at the LO (Born elastic scatter- 
ing) level, so that in each collision the transverse momenta of the two produced 
partons are equal, p±i = p±2 = P- This does not mean that these transverse 
momenta will be balanced in the rapidity window under consideration, so the 
event space for transverse energy deposition can be summarized by 

{0 if no particle gets into the gap 
p if one particle gets into the gap (3) 
2p if two particles get into the gap 

When considering the transverse energy production into a given rapidity window 
in pp collisions only the second and third possibilities are relevant. To quantify 
the computation of the contribution corresponding to cases 2 and 3 in (^ it is 
convenient to introduce the integral operators 

^1 = J dy,dy2i9{y^)+9{y2)-29{y,)9iy2)) *{...) (4) 
S2 = J dy^dy2{9{y^)9{y2)) *{...) , (5) 

where ^(^1,2) = 9{y„iin < 1/1,2 < ymax) Applying these operators to the differential 
cross section Eq. (|1|) we get a decomposition of the transverse energy production 
cross section in a given rapidity window into the separate one-jet and two-jet 
contributions (second and third entries in the event list Eq. (^)) 

dai da2 

+ 7^' (6) 



where 




(7) 



^In Eq. (H) and below pi = |p_Li 
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and 

da-?. 




(8) 

p=2E^ 



On the event-by-event basis these contributions correspond to completely dis- 
tinct possibilities of having the azimuthally balanced symmetric or unbalanced 
asymmetric transverse energy flow in the rapidity window under consideration. 

In Fig. (1) and Fig. (2) we plot the transverse energy production cross sections 
Eq. (0) and Eq. (g) for RHIC and LHC energies ^/s = 200 GeV and = 5500 
GeV, where for LHC we have chosen the energy to be available for protons in 
lead beams and the MRSG structure functions |^ were used. 

In Fig.(^ we also present the differential cross sections of transverse energy 
production into full rapidity intervals available at RHIC and LHC energies which 
will be used in the next section to normalize the differential cross sections of 
transverse energy production in pp collisions. 

The information contained in Fig. (1) and (2) is summarized in Table 1, where 
we show the parameters for the fits for the one-jet and two-jet spectra Eq. (|^) 
and Eq. (g having the functional form c{E±/l GeV)~". 



v^,GeV 


«i 


a2 


c, mb/GeV 


1 jet 


2 jets 


1 jet 


2 jets 


200 


5.09 


4.53 


173 


77 


5500 


4.25 


3.93 


3099 


819 



Table 1 



We see that at RHIC energies the one-jet asymmetric contribution dominates 
at small transverse energies, and at E± ~ 4.5 GeV the two-jet symmetric con- 
tribution takes over. At LHC energies the asymmetric contribution is clearly 
dominant in the whole minijet transverse energy range. 

3 Azimuthal asymmetry of minijet transverse 
energy flow in nuclear collisions 



In this section we turn to the analysis of the angular asymmetry of minijet 
produced transverse energy flow in nuclear collisions induced by the fundamental 
asymmetry in the pp collisions described in the previous section. The translation 
of the features characterizing the particle production in pp collisions to those 
characterizing the nuclear ones is possible, e.g., in a geometrical model, where 
the nucleus-nucleus collision is considered as a superposition of the proton-proton 
ones, see e.g. 0, 0. At each impact parameter b, where 6 is a distance in 
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Figure 1 One and two jets contributions to transverse energy production in pp 
collisions in unit central rapidity window at RHIC energy ^/s — 200 GeV. 




E,, GeV 



Figure 2 One and two jets contributions to transverse energy production in pp 
collisions in unit central rapidity window at LHC energy ^/s = 5500 GeV. 
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Figure 3 Transverse energy production in pp collisions in the whole rapidity 
interval at RHIC (v^ = 200 GeV) and LHC {y/s = 5.5 TeV) energies. 

the transverse plane between the centers of the colliding nuclei, the nucleus- 
nucleus collision is described as a Poissonian superposition of the nucleon-nucleon 
collisions such that a probability of n pp collisions is given by 

wnib) = ^NUb)e~'^--^'\ (9) 
n! 

where NAsib) is the average number of pp collisions in the nucleus-nucleus one, 
which is thus described as a specific superposition of multiple independent pp 
collisions occurring with the weight given by Eq. @. The elementary pp collisions 
occur between some nucleon belonging to the nucleus A located at the transverse 
distance bi from its center with the probability given by the probability density 
PA{bi) with the normalization 

J d%ipAib,) = 1. (10) 

and the nucleon from nucleus B located at the transverse distance 62 from its cen- 
ter with the probability given by the probability density ps(&2) with the collision 
probability P (b — bi + 62) • Thus the average number of pp collisions character- 
izing the basic Poissonian process Eq. (j^) is given by 

NAB{b) = AB j d%d%2 P{b-h + 62) PAibi) PB{b2), (11) 
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where P{b) is a probability of an inelastic collision of two nucleons initially sepa- 
rated by the transverse distance 6 — 61 +62- The physical meaning of the collision 
probability P{b) depends on the underlying physical mechanism responsible for 
inelastic transverse energy production in the binary nucleon-nucleon collisions. 
Our discussion is confined to minijets as providing such a source, so in our case 
P{b) is a probability of minijet producing inelastic nucleon-nucleon collision oc- 
curring at fixed impact parameter b. Let us stress that the differential probability 
of minijet-induced transverse energy production depends on the rapidity window 
under consideration. The usual assumption about the impact parameter depen- 
dence of the probability of nucleon-nucleon collisions P{b) is that the collisions 
are local in the impact parameter plane, i.e. 

m = <;™^'^*(|Ay|<yo)5(2H&)- (12) 

Let us now discuss in some details the normalization of the Poisson process Eq. (P) 
provided by the overall minijet production cross section into the rapidity window 
\^y\ ^ Voy cr™"*-''^*(| Ayl < Hq). The overall minijet contribution to transverse 
energy production cross section is given by the integral over E± of the differen- 
tial cross section Eq. (H). Because of the singular behavior of the perturbative 
transverse energy production cross section at small E± the very definition of 
the overall contribution of minijet mechanism to transverse energy production 
requires introducing a cutoff at small transverse energies 

Let us note that for any rapidity interval 

/ dE^^iAy) < a'-'(A|/), (14) 

where cr*"'^'(A?/) is an (experimental) inelastic cross section in a given rapidity 
window Ay. This shows that the cutoff Eq is physically a function of the consid- 
ered rapidity interval. 

Another important issue related to the choice of this cutoff is the possible 
contribution to the overall inelastic cross section of other mechanisms of trans- 
verse energy production, e.g. of the soft particle production due to the decay 
of hadronic strings. The restriction Eq. (|14D clearly refers only to the part of 
inelastic cross section corresponding to hard inelasticity, i.e. transverse energy 
production through semihard processes. The other part of the inelastic cross 
section corresponds to soft mechanisms of transverse energy production which do 
not involve large momentum transfers. It is important to note, that the char- 
acteristic time scale of semihard transverse energy production are smaller than 
that for the soft nonperturbative mechanism. At early stages of the collision hard 
parton skeleton is formed, which is then dressed by soft particle production due 
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to strings stretching in between the partons originating from primordial processes 
characterized by large momentum transfer. This shows, in particular, that the 
soft processes do not have, generally speaking, an independent share of the overall 
inelasticity, so the naive additivity 

da fi(jminijet ^^soft 

+ -nr- (15) 



dE± dE± dEi_ 

is, in general, not valid. It could happen, e.g., that with growing collision energy 
the yield of events with hard initial inelasticity would be dominant or even cover 
the whole event space (here we refer to non-diffractive contribution). In such 
extreme scenario the only function of soft mechanism is stretching the strings 
between the hard initial partons. Here it is important to recall that the cross 
section of transverse energy production Eq. (H), as computed in perturbative 
QCD, is a so-called infrared safe quantity and is thus entirely determined by its 
early quark-gluon stage and does not depend on the late stages of the process, 
including string formation between the separating partons. 

Let us stress once again, that in the present study we confine our consider- 
ation to analyzing the angular pattern of the primordial transverse energy flow 
generated at the early stages of the collisions by the semihard degrees of freedom 
(minijets). The analysis of the effects related to the subsequent redistribution of 
primordial transverse energy by soft interactions at larger times will be discussed 
in the future publications |]32| . 

The yield of perturbative contribution as a function of CMS energy is a cru- 
cial characteristics of the inelastic cross section. Unfortunately very little can 
currently be said about its magnitude, resulting in the uncertainty in fixing the 
cutoff for the perturbative contribution. 

In view of this we shall fix the cutoff value Eq at given collision energy as 
follows. To explore the possible "window of opportunities" for the hard minijet 
contribution as determined by the yield of independent soft particle production we 
will discuss two model scenarios. Namely, when considering the inelastic particle 
production at all rapidities we shall either assume the constant soft contribution 
a'^°^*{pp) = 32 mb, the inelastic cross section of pp scattering at intermediate 
energies, universally present at all CMS energies (mixed scenario) or assume 
that in all collisions the transverse energy is produced via the early perturbative 
minijet stage, i.e. put a'^°^^{pp) = (hard scenario). The cutoff Eq is thus 
determined from |^ 

cio"™"*-'^* ft^rrf f CTexp 5 soft contribution 

^inel 
exp 

(16) 



Eo ^ dE\ 1 0-*'^^' - 32 mb, soft contribution 32 mb. 



^The inelastic cross section is computed using the parameterization <Jinei{s) — ctq • 
(s/so)° °®''^ • (0.96 - 0.03 • log(s/so)), where sq = 1 GeV, ctq = 21.4 mb, which gives a good 
description of the existing experimental data |Q , see also a compilation in [|| 
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where the transverse spectrum in Eq. (|TBp refers to the full kinematic interval (see 
Fig. (^). Let us stress that the differential cross section for transverse energy 
production that we use in Eq. (|16D is the result of the lowest order calculation from 
the previous section, and the higher order effects that can phenomenologically 
be included within the geometrical unitarization scheme, see e.g. 0], are not 
included. Numerical values of the cutoff found by integrating the spectra shown 
in Fig. are given in Table 2 



^/S, GeV 


cr""^*, mb 




Eq, GeV 


Pi 


(T™"^^'"*, mb 


^PbPb , mb 


200 





41.8 


2.4 


0.54 


2.4 


5336 


32 


9.8 


3.5 


0.48 


0.54 


4102 


5500 





66.3 


6.9 


0.65 


2.8 


5443 


32 


34.3 


8.4 


0.64 


1.5 


4970 



Table 2 



In the fifth column we show the overall probability of asymmetric one-jet 
contribution pi{Eo) calculated using the differential transverse energy production 
spectra Eq. (D, (0) and (1): 



f oo 



ME,)^[l,E,^] rdE,4^]. (17) 



Ieo dE±l/ \Jeo dE 



As already mentioned, although the differential spectra describe the transverse 
energy production into some given rapidity window, in what follows the value of 



the cutoff Eq will be determined from Eq. (|Tj) considered for full rapidity window 
kinematically available for inelastic energy production at given collision energy. 
For more accurate determination of the cutoff Eq one would need experimental 
data on inelastic cross sections in, e.g., central rapidity window. Different quan- 
tities have different sensitivity with respect to the choice of the cutoff and, in 
particular, the dependence of Pi{Eq) in Eq. (|T^) on Eq is quite weak. In the last 
column in Table 2 we show the cross section of producing at least one minijet in 

minijt 
PbPb 



1 11 1 IT • minijet 

lead-lead collisions cr 



where WQ{b) is a probability of having no minijet producing nucleon - nucleon 
collisions, cf. Eq. (|^). 

The transverse energy production in nucleus-nucleus collisions is then de- 
scribed by the convolution of the distribution over the number of pp collisions 
obtained from Eq. @ at given impact parameter with the distributions charac- 
terizing the transverse energy production in pp collisions Eq. (|^ and Eq. (^. 

In practice this convolution was realized by a Monte Carlo procedure, where 
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• a large number (10^) nucleus- nucleus collisions were generated with the 
number of pp collisions distributed according to Eq. (P), 

• the weight of one-jet asymmetric (two-jet symmetric) pp collisions is equal 
to the probability pi (respectively 1 — pi) with pi taken from Table 2. 
More explicitly, this corresponds to a binomial distribution in the number 
of asymmetric collisions A^^^ 

^(iV,) = C^"pf (l-pi)« (19) 



• the weight for E± itself was in turn determined by Eq. (0) and Eq. (j|) for 
asymmetric and symmetric contributions correspondingly, 

• the azimuthal orientation of jet(s) was determined at random corresponding 
to a flat distribution in the azimuthal angle. For two-jet events the jets are 
going into opposite directions, so that their azimuth differs by vr. 

Let us now turn to the quantitative analysis of the event-by-event asymmetry 
of the minijet generated transverse energy flow. Our analysis will be made using 
a (normalized) difference between the transverse energy flow into the oppositely 
azimuthally oriented sectors with a specified angular opening 6ip each and rapidity 
window \y\ < 0.5. Let us note that this quantity has an important advantage of 
allowing for the future next-to-leading order analysis. For convenience one can 
think of the directions of these cones as being "up" and "down" corresponding 
to some specific choice of the orientation of the system of coordinates in the 
transverse plane. All our results are, of course, insensitive to the particular 
choice. Let us denote the transverse energy going into the "upper" and "lower" 
cones in a given event by E^{6(p) and Ef^{6ip) correspondingly. The magnitude 
of the asymmetry in transverse energy production can then be quantified by 
introducing a variable 

5E = E^{5ip)-E^{5ip) (20) 

Using the distribution over the number of asymmetric collisions Eq. ([T9|) and 
taking into account that the event space of asymmetric collisions is further sub- 
divided into two sets corresponding to nonzero energy going into the upper and 
lower cone E^ and we can calculate the quadratic mean of 6 E {6(f)) for the 
considered azimuthal openings 6(f) = 7r/2" {n = 0, 1, 2) 0: 

^/W) = ^eJp,N^^, (21) 
* 2 2 y ai — 6 

where ai is given in Table 1. The quadratic mean 6E{6(f)) in Eq. (|2ll) characterizes 
the magnitude of the disbalance in the minijet-generated transverse energy flow. 

■^It is easy to see that the distribution of 5E is a so-called multi-Poisson one 
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Note that 6E is essentially sensitive to the overall magnitude of the semihard 
(minijet generated) transverse energy flow. In Eq. ( pT]) this is clearly seen from 



{SE"^) oc Eq N. Numerical values for ^J {6E'^) are presented in Table 4. From now 
on we conflne our discussion to central PbPb collisions. 

In Figs.(^)-(|^) we show the probability distribution for SE in central PbPb 
collisions for two values of CMS energy 200 GeV and 5.5 TeV and two choices 
for Eq corresponding to mixed and hard scenarios. The angular apertures were 
chosen to be vr, 7r/2 and tt/A. From these flgures we see that for all types of 
collisions (except for Fig. (^) there appear peaks in probability distribution at 
values 6E = hEq. This is a reflection of a sharp cutoff adopted in the model 
and rapid decrease of minijet cross section with increasing E^. In the majority 
of cases this effect is seen only for small values of angular opening. The crucial 
parameter related to the peaks appearance is in fact N . The smaller is N the 
more evident are peaks. One could expect that hadronization and soft processes 
accompanying minijet production smoothed away these peaks. Curves that are 
initially smooth undergoes Gaussian law with dispertion {SE"^). Therefore, we 
can imagine appearance of curves with peaks after smoothing as Gaussian with 
dispersions given by Eq. (^l]) . For Fig. (|^) the result of such smoothing is shown 
in the inserted plot. 

Another useful quantity is a normalized asymmetry which is, on the contrary, 
insensitive to the absolute magnitude of transverse energy flow: 

E^{6ip) + Ei{6ip) 

where r G [—1,1]. In particular, the normalized asymmetry r simplifles compar- 
ing the asymmetries at different CMS energies. The values of r{6(f) in different 
collisions are characterized by the normalized probability distribution 



1 da 

P{r)\s^ = --r 
a ar 



(23) 

Sip 



To evaluate p{r) we use a Monte-Carlo simulation of the process of nuclear 
scattering as described above for the generated ensemble of 10'' PbPb collisions at 
RHIC and at LHC energies. We have calculated the asymmetry distributions p{r) 
for the central (zero impact parameter 6 = 0) collisions and cone apertures tt, 7r/2 
and 7r/4. The resulting probability distributions are illustrated in Figs. (|D-(|ll|) 
for mixed (ct"°^* = 32 mb) and hard (a'°^* = 0) scenarios at RHIC and LHC 
energies: 

Let us note that, in particular in the cases where the number of contributing 
collisions is not large (RHIC), one encounters "singular" conflgurations, for which 
r = —1,1 and r = corresponding to absolutely asymmetric and absolutely 
symmetric events in PbPb collisions. These are the events in which only one one- 
jet event contributes to the given aperture during the collision (r = —1, 1) or one 
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Figure 4 Probability distribution for the azimuthal asymmetry 5E_[_ in unit cen- 
tral rapidity window at RHIC energy ^/s = 200 GeV for central PbPb collisions, 
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Figure 5 Probability distribution for the azimuthal asymmetry 5E_[_ in unit cen- 
tral rapidity window at RHIC energy ^/s — 200 GeV for central PbPb collisions, 
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Figure 6 Probability distribution for the azimuthal asymmetry 6E± in unit cen- 
tral rapidity window at LHC energy = 5.5 TeV for central PbPb collisions, 
crsoft ^ 
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Figure 7 Probability distribution for the azimuthal asymmetry 6E± in unit cen- 
tral rapidity window at LHC energy -^/s — 5.5 TeV for central PbPb collisions, 
a'^f* = 32 mb. 
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Figure 8 Probability distribution of the normalized azimuthal asymmetry p{r) 
in unit central rapidity window at RHIC energy ^/s = 200 GeV for central PbPb 
collisions, cr*"-^* = 




Figure 9 Probability distribution of the normalized azimuthal asymmetry p(r) 
in unit central rapidity window at RHIC energy ^/s — 200 GeV for central PbPb 
collisions, o"*"-^* = 32 mb 



16 
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Figure 10 Probability distribution of the normalized azimuthal asymmetry p(r) 
in unit central rapidity window at LHC energy y/s = 5.5 TeV for central PbPb 
collisions, cr*"-^* = 




Figure 11 Probability distribution of the normalized azimuthal asymmetry p{r) 
in unit central rapidity window at LHC energy ^/s = 5.5 TeV for central PbPb 
collisions, a"^"-^* = 32 mb. 
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two-jet event contributes to r = 0. Their probabilistic weight can be described 
hj 6 - functional contribution to p{r) at the "singular" points. Their yield in the 
minijet event ensemble is given in Table 3 for mixed = 32 mb) scenario at 

RHIC energy (other values are negligible). 



6{p 


r = —1 


r = 


r = 1 


tt/2 
tt/A 


1.3-10-3 
2.5 - 10-2 


1.6 - 10-2 
1.2 - 10-^ 


1.3 - 10-3 
2.5 - 10-2 



Table 3 



In Fig. these contributions would correspond to infinitely narrow peaks 
and are thus not shown. 

The angular pattern of the transverse energy production as characterized 
by the considered energy-energy azimuthal correlation probability distribution is 
conveniently described by the lowest moments of p{r). In Table 4 we present, 
together with the numerical values of the quadratic mean ^ {SE'^), cf. Eqs. 
and (pTj), the values of the standard deviation a defined as 

a2 = j dr {r - ff p{r) , (24) 

where in our case (r) = 0: 



(jSoft 

mb 


GeV 


NpbPb 


GeV 




a 


6(j) = Ti 


50 = 7r/2 


50 = 7r/4 





200 


75.7 


21 


0.084 


0.088 


0.124 


0.178 


5500 


87.4 


84 


0.086 


0.094 


0.133 


0.189 


32 


200 


17.1 


14 


0.168 


0.177 


0.259 


0.387 


5500 


47.1 


74 


0.117 


0.127 


0.180 


0.257 



Table 4 



Note that all the data presented in Table 4 include contributions from the 
singular points r = —1, 0, 1. 

Let us now turn to the analysis of the results presented in Figs. (4)- (11) and 
Tables 3-4. The main goal is to understand the dependence of the angular pattern 
of the transverse energy flow on the basic parameters such as the infrared cutoff 
Eq, total number of minijet-generating collisions iV, the yield of asymmetric pp 
collisions pi and CMS energy a/s. 

In the forth column in Table 4 we show the quadratic mean {SE"^) for the 
azimuthal opening 50 = vr. The results agree with Eq. (pl|), so that the average 
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disbalance in the transverse energy is indeed essentially determined by Eq and 
N. 

To understand the results for normalized asymmetry p{r) it is helpful to con- 
sider a simplified model, in which the elementary pp collisions can produce only 
some given amount of transverse energy, 

/ \PP 

' ^ ^ = a'^-'^(v/^)5(E^-Eo(v^)), (25) 



so that all transverse energy is assumed to be produced exactly at the cutoff Eo. 
Note that except for ascribing the energy production to elementary pp collisions 
this model is very similar to the expected pattern of transverse energy production 
in the quasiclassical approach based on McLerran-Venugopalan model, cf. 
Then, for considered azimuthal apertures 6(j) = 7^/2^ [n = 0, 1, 2) we get for the 
standard deviation a defined in Eq. 




-[^^=^)^\l:^:-^r{^ + 0[-^)) (26) 

This shows that the width of the distribution p{r ) is determined by the ratio of the 
relative yield of asymmetric collisions pi to the average number of collisions. In 
Table 4 we compare the predictions of this simple model to the values of standard 
deviation a computed using the differential spectra plotted in Figs (1), (2) (to 
save space, only the results for 50 = vr are given in column 5) and observe only 
a 10 % difference. This shows that the results obtained using the continuous 
spectra in Figs. (1), (2) are essentially determined by the contribution at the 
cutoff energy Eq. 

From Fig. (^ we see that for small number of asymmetric collisions the shape 
of p(r) has peculiar sharp peaks at certain values of r. The origin for this is in fact 
the growth of the differential cross section for transverse energy production at 
small E± in pp collisions, cf. Figs. 1-2. Indeed, let us, for simplicity, assume that 
each pp collision in the restricted minijet ensemble can produce the transverse 
energy exactly at the cutoff E± = Eq only, cf. Eq. (p5[) . In this case in addition 
to the "true" singular points r = —1, 0, 1 we will have " semisingular" ones so that 
for a particular event containing n minijets, with ni| being the number of "up- 
coming" one-jet events, being the number of "down-coming" one-jet events, 
and n2 = n — ni^ — being the number of two-jet events, the following exact 
relation does hold: 

r = (27) 

n 

Thus the values of r belong to a set of rational numbers in the interval [—1,1], 
which we call "semisingular". Of course, the most spectacular "semisingular" 



The details of this calculation can be found in Appendix A 
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points are those with small numerators and denominators both due to the higher 
frequency of events having small number of minijets and to smaller distribution 
width (deviation from E± = hEq) for events with small number of asymmetric 
collisions. 

Let us note that the appearance of singular points —1, 0, 1 is a consequence of 
calculating the cross sections for transverse energy production in the elementary 
hard block in the lowest order in perturbation theory. In the next-to-leading 
order, where the transverse energy can be shared between three (mini)jets these 
singular points will become milder singularities of p{r) at r = —1, 0, 1. This shows 
that the calculation of the true shape of p{r) near the singular points requires, 
as usual, resumming the perturbative contributions to all orders. 

Physically, within the scheme adopted in this paper, the number of semihard 
collisions depends on the sharing of inelasticity between soft and hard mechanisms 
of transverse energy production. In the mixed scenario we assumed that 32 mb 
of the inelastic cross section of pp collisions corresponds to the soft production 
mechanism, while the rest of inelastic cross section is due to semihard production. 
In the hard scenario it is assumed that the semihard transverse energy production 
saturates all available inelasticity. It is to be expected that the distributions 
characterizing azimuthal asymmetry of transverse energy production p{r) defined 
in Eq. ( p3D will be wider in the mixed scenario than in the hard one. From Table 5 
we see that this is indeed the case. At RHIC energies the standard deviation for 
mixed scenario is bigger than that in hard one in the ratio acinic /^'rh'ic — 
2.0 — 2.1. At LHC energies the effect is less pronounced, here a^jfcf/^LHC — 
1.3 — 1.4. We see that with growing CMS energy the angular pattern of the 
minijet-generated transverse energy flow is becoming less sensitive to the relative 
weight of perturbative and nonperturbative contributions to the inelastic cross 
section. 

The dependence of the standard deviation a on the aperture, remains es- 
sentially the same for both CMS energies and values of the impact parameter 
considered and is inversely proportional to the angular opening: 

an/2n ^ 22 (28) 

which is consistent with the prediction of a simple model of transverse energy 
production Eq. (^BD and corresponds to purely statistical change in the standard 
deviation, where shrinking the angular aperture by a factor of 2 enlarges the 
standard deviation by a factor of ^/2. 



4 Conclusions 

The main results of our analysis can be formulated as follows. 

We first discussed a basic asymmetry in the minijet transverse energy pro- 
duction in a restricted rapidity window in pp collisions due to different proba- 
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bilities of having a "symmetric" two-jet or "asymmetric" one-jet contribution in 
the rapidity interval under consideration. The cross sections for symmetric and 
asymmetric contributions in pp coUisions for RHIC and LHC energy show that 
while at RHIC energy the weight of both configurations is approximately equal, 
at LHC energy the asymmetric contribution is clearly dominant. 

We further considered a geometrical model for nuclear collisions in which they 
are described as an incoherent superposition of nucleon-nucleon ones. We dis- 
cussed two possible partitions of the inelastic cross section in terms of soft and 
semihard contributions and analyzed the angular pattern of minijet-generated 
energy flow for central and peripheral nuclear collisions at RHIC and LHC en- 
ergies. Specifically we considered probability distributions for transverse energy- 
transverse energy correlations in the oppositely azimuthally oriented cones with 
varying aperture. We show that the resulting distributions are essentially sen- 
sitive to the number of semihard collisions, which is in turn dependent on the 
above-mentioned partition of the inelastic cross section into contributions of dif- 
ferent types, and on the (related) choice of the infrared cutoff. We also show that 
the results are very close to the predictions of a simple model, in which all the 
transverse energy is produced directly at the infrared cutoff. 

The approach developed in this paper could be further generalized to the 
analysis of the minijet generated background oriented flow (for the definition 
of the oriented flow and comprehensive discussion see e.g. |^). In particular, 
as the importance of minijet contribution is expected to grow with energy the 
presence of the background oriented flow of purely fluctuational origin could 
increasingly influence corresponding hadronic observables. 

Another crucial issue is the dynamical evolution of the primordial partonic 
inhomogenities in the course of parton-hadron conversion. In a recent study |^ 



it was shown, that the seed inhomogenity in the initial condition of the elliptic 
flow type for the hadronic RQMD code survives the freeze-out and is visible in 
final azimuthal distributions. This question is surely most important and will be 
discussed in the forthcoming publication p2|. 



Acknowledgements 



We are grateful to K. Kajantie, G. Zinovjev, J. Schukraft, K. Eskola, I.M Dre- 
min, I.V. Andreev, J. -P. Blaizot, J.-Y. Ollitrault, L. McLerran and P. Jacobs for 
useful discussions. 

A.L. is grateful for kind hospitality and support at CERN Theory Division, 
where this work was started and at Service de Physique Theorique de Saclay, 
where the major part of it was done. 

We are grateful to the referee of the first version of this paper for constructive 
remarks and suggestions. 

This work was also partially supported by Russian Fund for Basic Research, 
Grant 96-02-16347. 



21 



The work of D.O. was partially supported by the INTAS grant 96-0457 within 
the ICFPM program. 



References 

[1] S.A. Baas, M. Gyulassy, H. Stoecker and W. Greiner, "Signatures of Quark- 
Gluon Plasma Formation in High Energy Heavy Ion Collisions", | |hep-| 
I ph/981028ll 



[2] E.M. Levin and M.G. Ryskin, Phys. Repts. 189 (1990), 267 
[3] X.-N. Wang, Phys. Repts 280 (1997), 287 

[4] K.J. Eskola, Comments. Nucl. Part. Phys. 22 (1998), 185 [^iucl-th/9705027| ]; 
K.J. Eskola, "Minijets in ultrarelativistic heavy ion collisions at RHIC and 
LHC", Proc. RHIC Summer Study '96, Theory Workshop on Relativistic 
Heavy ion Collisions, July 8-19, 1996, Brookhaven National Laboratory, NY, 
USA, eds. D.E. Kahana and Y. Pang [|hep-ph/9610365| ; 

[5] K. Kajantie, P.V. Landshoff and J. Lindfors, Phys. Rev. Lett. 59 (1987), 
2527 

[6] K.J. Eskola, K. Kajantie and J. Lindfors, Nucl. Phys. B323 (1989), 37 

[7] K.J. Eskola, K. Kajantie and V. Ruuskanen, Eur. Phys. J. CI (1998), 627 

[8] K.J. Eskola and K. Kajantie, Zett. Phys. C75 (1997), 515 

[9] K.J. Eskola, K. Kajantie, P.V. Ruuskanen and K. Tuominen, "Scaling of 
Transverse Energies and Multiplicities with Atomic Number and Energy in 
Ultrarelativistic Nuclear Collisions" | hep-ph / 9909456] 

[10] T. Sjostrand and M. van der Zijl, Phys. Rev. D36 (1987), 2019 

[11] J. P. Blaizot and A.H. Mueller, Nucl. Phys. B289 (1987), 847 

[12] L. McLerran and R. Venugopalan, Phys. Rev. D52 (1994), 2233; D49 (1995), 
2225 

[13] A. Kovner, L. McLerran and H. Weigert, Phys. Rev. D52 (1995), 3809; ibid. 
D52 (1995), 6231 

[14] Yu. V.Kovchegov and D.H. Rischke, Phys. Rev. C56 (1997), 1084 

[15] S.G. Matinyan, B. Mueller and D.H. Rischke, Phys. Rev. C56 (1997), 1927; 
ibid. C57 (1998), 2197 



22 



[16] M. Gyulassy and L. McLerran, Phys. Rev. C56 (1997), 2219 



[17] Yu. Kovchegov and A. Mueller, Nucl. Phys. B329 (1998), 451 
[18 



A. Krasnitz and R. Venugopalan, "Real time simulations of high en- 
ergy nuclear collisions", | hep-ph /9808332 1 ; " Non-perturbative computation 
of gluon minijet production in nuclear collisions at very high energies" 



hep-ph/ 9809433|| ; Making glue in high energy nuclear collisions", [[hep- 



ph/9905319|| ; "The initial energy density of gluons produced in very high 
energy nuclear collisions" , |Piep-ph/ 9909203| 



[19] S.A. Bass, B. Mueller, W. Poeschl, " Lattice Gauge Description of Colliding 
Nuclei", [nucl-th/980811] 

[20] W. Poeschl, B. Mueller, "Gluon Field Dynamics in Ultra-Relativistic Heavy- 
Ion Collisions: evolution on a Gauge Lattice in 3-1-1 Dimensions", [|nucl-| 
th/9808031 |; "Gauge Lattice Simulation of the soft QGP Dynamics in Ultra- 
Relativistic Heavy- Ion Collisions", ||hep-ph /981 1441|] ; "Real-Time Evolution 
of Soft Gluon Field Dynamics in Ultra-Relativistic Heavy- Ion Collisions", 
|nucl-th/ 9903050 1 



[21] 

[22] 
[23] 
[24] 



A.H. Mueller, "Toward Equilibration in the Early Stages After a High En- 
ergy Heavy Ion Collision" , | [hep-ph / 9906322|] ; " The Boltzmann equation for 



gluons at early times after a heavy ion collision" , ||hep-ph/990938 



L. McLerran, "Three Lectures on the Physics of Small x and High Gluon 
Density", [[hep-ph/9903536[| 



A.H. Mueller, "Small-x Physics, High Parton Densities and Parton Satura- 
tion in QCD", |[hep-ph/ 99 112891 



R. Stock, "Analysis of Single Events in Ultrarelativistic Nuclear Collisions" 
Proc. Pre-Conference Workshop "Quark Matter 95", [UCRL-TD-121571] 



[25] E. Shuryak and CM. Hung, Phys. Rev. C57 (1997), 397 
[26] 



[27] 
[28] 



M. Gazdzicki, A. Leonidov and G. Roland, Eur. Phys. Journ. C6 (1999), 
365 

M. Stephanov, K. Rajagopal, E. Shuryak, " Event-by- Event Fluctuations in 
Heavy Ion Collisions and the QCD Critical Point" , |[hep-ph/ 9903292 



A. Leonidov and D. Ostrovsky, "Minijet Transverse Energy Production in 
Next-to-Leading Order in Hadron and Nuclear Collisions" | hep-ph/9811417 
to appear in Eur. Phys. Journ C 



M. Gyulassy, D. Rischke and B. Zhang, Nucl. Phys. A613 (1997), 397 



23 



[30] X.-N. Wang and M. Guylassy, Phys. Rev. D44 (1991), 3501; D45 (1992), 
844; Comp. Phys. Commun. 83 (1994), 307 



[31] A.D. Martin, R.G. Roberts and W.J. Stirling, Phys. Lett. B354 (1995), 155 

[32] A. Leonidov and D. Ostrovsky, in preparation 

[33] UA4 Collab., M.Bozzo et al., Phys. Lett. B147 (1984), 392; 
UA5 Collab., G.J.Alner et al, Z. Phys C32 (1986), 153; 
Fermilab E710 Collab., N.Amos et al, Phys. Rev. Lett. 63 (1989), 2784; 
R.M.Baltrusaitis et al., Phys. Rev. Lett. 52, 1380 (1984) 

[34] A. Leonidov, J.-Y. Ollitrault and D. Ostrovsky, in preparation 



[35] J.-Y. Ollitrault, "Flow Systematics from SIS to SPS energies" ||nucl- 



I ex/ 9802005 [ 

[36] F. Wang and H. Sorge, "Remnants of Initial Anisotropic High Energy Den- 
sity Domains in Nucleus-Nucleus Collisions" ||nucl-th/9811006[ 



Appendix A (r^) calculation 

In this Appendix we present a derivation of the formula for the standard deviation 
a Eq. dH). 

Let n be the number of minijet producing hadron-hadron collision in a given 
the nucleus- nucleus one characterized by the Poissonian distribution Eq. (P). Let 
us further denote by ni the number of those minijet producing hadron collisions 
in which only one minijet hits the rapidity window under consideration (asym- 
metric contribution), and by rii^ and are the numbers of such single minijets 
propagating into the upper and lower of the two oppositely oriented cones corre- 
spondingly. In what follows we shall use a simplified model of transverse energy 
production in which it is produced strictly at the cutoff Eq (cf. Eq. (pSf). The 
averaging over event ensemle has to be done in the sequence opposite to the 
one adopted in the Monte Carlo procedure. First, we average over rii^ at fixed 
ni = + ni|: 

ni|=0 

Next, we average over rii at fixed n according to the binomial probability distri- 
bution Eq. (0): 
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Finally, we have to average over the Poissonian distribution Eq. yielding 

(A 

Where F{1, 1; 2, 2; A^) is a generalized hypergeometric function. 
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